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ilBSTHACT 


In this dissertation^algoritiims ensuring optimal 
results for following three problems concerning locations 
on networks are given. 

(1) Finding all m-vertex median of a graph, 

(2) Finding CBlmodified m-vert ex- median of a graph. 

( 3 ) Finding the minimum number of emergency service 
facilities required for a network of destinations 
under certain assumptions. 

For the first two problems only heuristic methods or 
complete enumeration are available. lo find the optimal 
results for the above problems, the distinct approach adopted 
in this dissertation is to reduce the above problems to mini- 
mization of nonlinear pseudo-Boolean objective functions with- 
out any constraints, 

Basic algorithm is available to minimize the unconstrained 
nonlinear pseudo-Boolean functions. 

Computer codes are developed for the algorithms proposed . 
here and computational attributes of various algorithms are . 
compared. 
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CHiPTER I 


IRTRODirOTIOR 

> Jc- . 

Determining good locations for sources on a network 
has received a reasonable amount of attention in the last 
decade. Some of the in^jortant problems, regarding optimal 
locations of these sources are: 

(1) The problem of supplying n destinations from m sources 
has been attacked with a variety of assumptions and methods. 

If both sources (with given capacities) and destinations 
(with given demands) are at fixed locations, then this becomes 
a standard transportation problem. 

But here the task is to determine the locations of the 
sources so as to minimize the transportation cost, when the 
locations for the destinations are given with certain assum- 
ptions. If the destinations consist of fixed vertices on a 
network, but the sources may lie any where on the network-* 

edges and the destination demands are fixed, while the source 

( 2 ) 

capacities are unconstrained, it has been shown by Hakim i ^ 
that the problem resolves itself (by taking transportation 
costs proportional to distance) into finding the generalized, 
absolute median of the corresponding graph to the network in 
which both vertices and the edges have weights. The weight 
of a vertex representa a destination demand and the wei^t 
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of ^edge represents tiie shortest distance "between the 
corresponding vertices in the netvi/ork. Hakimi also demonstra- 
ted that there will exist such ah. n ^ . median that includes 
only vertices of the graph. Thus a solution to this problem 
will correspond to the case where both destinations and the 
sources lie on the vertices of a network, a situation like 
that investigated for fixed destination demands and the un- 
constrained source .capacities by Maranzana^®\ 

(2) While locating the sources, an important consideration 

is the maximum distance that separates a destination from its 

nearest source. In first problem, by imposing a limit on the 

a 

maximum distance that seperates^destinat ion from its nearest 

source, we will have to solve the first problem under n new 

constraints (one for each destination and hence in corres- 

» 

ponding graph G- caae for each vertex) to meet the specified 
distance standards. This is exactly what the second problem 

ClG ) 

is and has been named as ’Modified m-vertex median of a graph'\ ■ 

(3) In the third problem we have to find the minimum number 
of sources required* to meet the distance standards .or 
constraints (i.e. limit on the maximum distance that seperates 
a destination from its nearest source) of each of the desti- 
nations . As such this problem is especially applicable 

to the location of emergency service facilities (consider 
sources as service facilities here) such as fire stations , 

■^he assumptions made are same as in the previous problems. 
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although one may equally well apply it to the location of 
ordinary services, such as schools, libraries etc* 

.The location of fire stations might be approached 
according to the structure just described. The limit on 
response time/distance (as time may be taken proportional to 
distance) is imposed to ensure, that no more than a specified 
time period s^ will elapse before a response will occur to 
any call at jth destination. The desired solution to this 
problem locates the minimum number of fire stations that 
satisfies the response time requirements. 



OHiPTER II 


A SURVEY OP MLISIIK& ALGORITHMS 

In this chapter, a surrey of the existing heuristic 
methods is done for the first prohlem. The heuristic methods 
for the first prohlem are modified so that they can he used 
for the second prohlem. Also the optimal algorithms available 
for the third prohlem are mentioned, 

( 2 ) 

The m-Vertex Median ^ ■ 

Consider a graph G- with weighted vertices and edges* 

SBie iL<ytatioiiE used axe as follows: 


Y 

Yertex set of G having n vertices v^ , , . 

, V . 

> n . 


- Weight of the vertex v. 


A3 

- Weight of the edge(i, 3 )i, e. the.shortest 

distance 


between vertex v. and the vertex v^ • 

X cl 


L 

- Distance matrix (symmetric) with dj^^ *s as 

the 


elements. 



If in the graph G, all vertices have identical wei^ts, ’ 
a vertex median Vj^ will he that vertex for which the sum 
of the elements in the corresponding column of I) is minimized. 
That is 5 let 

n ^ 

d. = 2 d. . (for 3 = 1 

3 ^ 


i ♦ • 


n) 


( 2 , 1 ) 
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tlaen is the vertex median if and only if 

djj. = min (d-j , . . . .,d^) (2.2) 

If in the graph G, the vertices have unequal weights, 
then it is necessary to redefine the distance matrix, let 
2 be the nth order diagonal matrix with the vertex weights 
on the diagonal. Now the weighted distance matrix of G 
is defined by 



(2.3) 


Matrix H is no longer symmetric. Each h. . represents 

1 J 

the wei^ted distance associated with vertex v. if v. were 

J o. 


the unique source. The vertex median is defined as above in 

* 

Eqn, (2.1) and Eqn. (2,2) but over the matrix H. 


1 generalization of vertex median follows logically 
from Eqn. (2,l) and Eqn. (2,2). Let be a subset of Y 
containing exactly m vertices of G. Eor an n vertex graph 
there will be (^) possible subsets of cardinality m and we 
arbitrarily index them "by 7^ (p = 1,2, , , , (^))» each 

subset we may construct a submatrix of H by adjoining 
all rows of H for which the corresponding vertices are 
contained in 7^. is of order mxn and describes precisely 
the sources and the associated weighted distances for every 
destination if the set of sources is limited to vertices in 


* Henceforth whenever we use Eqn, , we mean either equality, 
inequality, or syst an of inequalities and/or equalities. 
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it is assumed that sources have no capacity 
constraints, each destination v^ will be served by that 
source v^^ in for which hj^^ is a minimum, i.e. 




1 m 


( 2 . 4 ) 


The total wei^ted distance hp for. the set of sources 
will be the sum of column minima of 



n 


3=1 


3 


( 2 . 5 ) 


where k refers to the source for which h. . is minimized, ■ 
An m vertex median of G- is now defined as same such that 


hp = min (h^ , h^ , 


so that hp <_ hp (p = 1 , 2, 



( 2 . 6 ) 


-vertex 

The m median- is not necessarily unique. Clearly it defines a 
A 

set of sources which is in some sense closest to all desti- 
nations. 

So far a completely satisfactory method for finding the 
m-vert ex median of a graph is not available. However, two 
heuristic approaches available till now are described below. 
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(s ) 

2,1*1 Maranzatia’ s method ^ 

The algorithm starts with an arbitrary selection of 
a subset 7^ containing m vertices and partitions the graph 

into m subsets such that to each vertex in 7 . one subset 

m ’ 

of these m subsets is associated. This is accomplished by 
associating each vertex in G- with its nearest vertex in 7^^. 
Next the center of gravity of each set in the partition is 
determined and the original vertices in 7^ are replaced by 
these vertices. The process is repeated until the vertices 
in 7jj^ do not change, A formal statement of algorithm follows: 

Algorithm 

Step 1: Arbitrarily select m vertices in 7 and assign these 

vertices to the variable array, 7 appearing in 

^i 

Step 2, 


Step 2; Associated with this array of m vertices, v. 






V determine a corresponding 
^m 


*0enter of gravity: 7ertex Vj^^ will be a center of gravity 

t 

of a subset 7 of 7 if 


Z 

V 


r 



< 2 i 

v.e7 



for all ■Wjj.67 
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partition of P_ . "by putting 

^1 ^2 


^yi = i I \ 


> k i V kf°^ yj 


If a vertex is equidistant from more than one v , a decision 

is required relative to placement of the vertex. lo break 

the tiey put the vertex in set associated with v having 

7 ± 

smallest among them. 

Step 3: Determine -tibecent re of gravity, n for each P.. , 

Step 4: If computation is stopped and 

the current values of v„ and ?_ constitute the 


desired solutioni otherwise set v^ 
to Step 2, 




and return 


In Step 3 if the center of gravity is not unique, 
then choose one w'hich has smallest subscript among them* 

(8 

The algorithm has been shown to be monot onically' ' 

non-increasing' with respect to the successive selection of 

'P according to Step 2 and successive selection of value 
Ti 

of 0„ for Py according to Step 


19 I 

2.1,2 One Optimal or Vertex Substitution Method 

eons id er the definition of the m vertex-median by 
Eqn. (2.5) and Eqn. (2.6). Por each possible subset of 
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vertices we may construct a sutmatrix of H ty adjoin- 
ing m rows corresponding to the vertices v. in . The vertex 
in with which any vertex v^ is associated is defined as 
that Vj^ such that 

This is the jth column minimum in The total wei^ted 
distance h^ for the pth subset will be the sum of these 
column minima. 

Suppose that we decide to replace one vertex Vj^, in 
p 

the subset Yjj^ by another v-j^. Several kinds of effect on the 
total weighted distance may occur. 

If h.. were not the jth column minima, then its 

J. J 

replacement by h^. mi^t have several outcomes depending 


rJ 

on whether ^ 



(2.8a) 

or 


^tj ^ ^3 

(2.8b) 

where h^j is jth column minima of 



In case of Bqn. (2,8a), the jth column contribution to h 
will be zero i. e. j '^bi = 0., 
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In case of Eqn. (2,8'b) 


cindl 


3 ti = 

D bi < 0 


If were the jtb. column minima, following will be observed; 

(2.9a) 
(2.9b) 
(2.9c) 

where h . is the wei^ted distance from vertex v. to that 

S3 3 

vertex v for which 
s 

h. j for e 

(2. 10) 

and Vg* ^ 

In otherwords, h^- is the second smallest nth column 
el ement in . In cas e of Eqn. ( 2. 9a ) , the 3 th column cont ri- 
but ion to h from the substitution of v-^^ for v^ is now incre- 
mented by 

3>^ bi = h^^ - h^ (2.11) 

3^:bi' < 0 " 




< ii- • 
DO - 13 


Or 


h. . < h, . < h . 

10 ” b3 - S3 


or 


1 < IIk-; 

10 ~ SO ■“ DO 


and 



In case of Eqn. (2<.9’b)5 the jth column contribution to h is 
incremented by 


jA-bi = h_. - h 

Ck-nd 

• A 

0 bi > 0 


( 2 . 12 ) 


In case of Eqn. (2»So), the jth column contribution to h is 
incremented by 


cxnd. 


.A 

J 

• A 

D 


bi = 

bi > 0 




(2*13) 


whether it is worth substituting vertex for depends : 
upon the net effect of the increments summed over all column 

l>i = S 3^ bi (2.14) 

Substituting v^ for reduces tota,l wei^ted distance only 
if 

^ bi ■ < .0 

These observations direct us to develop a . One optimal 
algorithm as follows. 
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Algorithm 

Step 1: Select some initial vertex subset V^jj.i'or convenience 
in expression let it contain vertices v^ , Vgs . ..-*5 

V j Vjj^. 

Step 2: Compute the total weighted distance h^^^ for the 
system. 


Step 
Step 4-: 

Step 5i 

and 


Select some vertex, v-j^, not in the subset Y^» 

"For each vertex v. in substitute v, and compute 
- 1 m b 

bi. 

Find that vertex in V_, such that 

A. m 

^bk < 0 (2.16) 

‘~^bk = min ^ bi f or (i = 1 , 2 , , , n') 


Step 6i If a vertex satisfying the Eqn. (2,16) can be found, 
substitute v-j^ for v^ in the subset Vj^, label the 

I ' 1 

new subset so formed as Y^^^ and compute h^j^. If no 
vertex satisfies the Eqn. (2,16), retain the subset 
Yj^ and proceed to Step 7, 

: 1 

Select another vertex not contained in Y^ or Y^ 
and not previously tried and repeat St eps 4 'through 


Step 7 : 
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Step 8 ! When all the vertices in the compliment of 

ft 

have been tried, define the resulting subset 
as new and repeat Steps 7 through 9. Call each 
such complete repetition a cycle. 

Step 9 ! When one complete cycle of Steps 3 through 8 

results in no reduction in hjterminate the procedure. 

fi 

The final is ail One optimal estimate of an 
m-vertez median of G-, 

At first sight this procedure appears laborious 
but it is much superior to Maranzana's method in many aspects 
and that will be quite evident in the Chapter Y of this 
dissertation. 


(■-in) ' 

2.2 The Modified m-Tertex Median ^ '' 

As mentioned in the Chapter I, this problem is 
nothing but the first problem with n new constraints, one 
for each vertex in the graph G-. 

If s . is the limit for the vertex v. i.e, the 

maximum weighted distance allowed that sepearates v^ 
from its nearest source vertex, -then set 
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will "be tile set of vertices within s. of v. and can provide 

J J 

acceptable service to the vertex v.. Jor n vertices there 

t] 

will he n such sets and no set will be empty as h.. is 

33. 

always zero while s . is a positive quantity, 

3 


Here one may notice that every subset containing m 
vertices of V may not be a feasible solution set to the 
problem as in the first problem. It may also happen, there 
exists no feasible solution set to the problem. 

Ary solution set will be feasible if and only if 

(2.18) 

for (3 =1,2,. , , , n) 


If it is not possible to get such a v\?hich satisfies 

I-O, 

the Eqn. (2,18), the problem will have no feasible solution. 

So far a satisfactory method for solving the problem 

is not available. In the following subsections, the existing 
heuristic methods to solve the first problem are modified to 


solve this new problem, 

2,2.1 Modified Mar anz ana's method: 

The starting solution taken for this algorithm 
should be feasible' too. This can be obtained by solving the 
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following system; 


i6 




3^1 > 1 


(2.19) 


for (j = 2, , . , , n) 


where 


n 

Z y = m 
i=1 ^ 


6 I Oj Ijand a 


= h I V. 






If (y^ 5 « • is any feasible solution then corresponding 

j Ai 

the feasible solution set V^n = ; y,- = 1 • < 

There are many methods to solve system (2.19). One 
elegant way is pseudo-Boolean method. 

The modified algorithm is given below; 


Algorithm 


Step 1; Get a feasible solution set as described above. 


let it contain m vertices v^ 


* 


* *’ ^yn 


St 


ep 2; Determine a corresponding partition of Y, JP^ ’ 

■^1 


J * ♦ f 


by putting 

% 


^yj = ^ ^y.lc and y. e 


step 3; Determine a centre of gravity, , for each P„ 

■^3 -3 

satisfying 1 ^'°n v^^ € Py also. 


*In the algorithm y^ is not bivalent variable but at subscript 
which can take values from 1 to n. 
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(where Vn = 0 ) 

Step 4: If 0„ = for all j? computation is stopped and 

the current values of v_ and P„ constitute the 

desired solution, otherwise set v_ = C and return 

■^3 ^3 

to St ep 2 , 


2.2,2 Modified One optimal or Vertex substitution method: 

Here also we have to start with a feasible solution 
as described above by solving system ■ (2,19). While looking for 
a substitute of a vertex we have to always check whether the 
resulting solution set is feasible or not and h^ for solution 
set is calculated as follows: 


and V 6 H . 

■ J 

Other computations are same as in One optimal method, 

2 , 3 binding Minimum Humber of Etoergency Service Facilities ^ 

In this problem we have to minimize the total number 

of sources required to meet the distance standards for, each 

( 2 ) 

of the destinations, Hakimi^ ' using Boolean functions, sought 
the minimum number of sources that covered all the destinations 


iU 


n 


° .1=1 '"'^3 


Where hj^.^ < , , Tj, 6 
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eve'^3 

witliin the specified maximum distance for that destination. 
The resulting method requires an enumeration of all feasible 
solutions and as the problem size grows, the effort of deter- 
mining the minimum number of sources can be expected to grow 
rapidly. 

The problem can be mathematically expressed as 

follows : 

n 

Min 2 y . 

1=1 ^ 

such that 

2 y. >1 (2.20) 

lea ^ “ 

for (o = 1 , 2, , . , , n) 


5 , a. have same definition as in the ‘second 

yd 


Where y. 0 , 1 
^ i 

problem. 

There are three other methods which seem most favoured , 
these are! 


( 10 ) 


(i) linear programming and cutting plane techniques | 

(ii) Reduction techniques 5 

(iii) Implicit enumeration. 

The fourth approach developed is described in 


Chapter 17 



GHiPTER III 


BASIC GONOEPIS OE PSEUDO-BOOIEilf PROGRAMING 


In this chapter "basic concepts of pseudo-Boolean 
programming are reviewed^^ 

Definition 3*1? 


Any real valued function with bivalent variables is 
called a pseudo-Boolean function ^- Thus f will be a pseudo- 
Boolean function if 


f : 



where R is the field of real numbers and ^2 ~ \ 

obivious that such a function can be written as a pol3m.omial 

linear in each variable. 


■Thus a function f(x^ , x^) can be written 

linear in ith variable as follows: 

f(x^> , . , , x^) = x^ g(x^ y .»»; ’ ^i+1 ^ ••• 


+ h(x^ j ^2* • • • 5 > ^i+i ^ 


(3.1) 


in equality of the form 
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where f is a pseudo-Bo clean fimction, is called a pseudo- 
Boolean equality , • 

If Eqn. (3*2) would have heen an inequality then it 
will be called pseudo-Bo clean inequality . 

Definition 3.5 ; 

i PS eudo-Bo clean programme is the problem of optimizing 
a pseudo-Boolcan function, whose variables can be either un- 
constrained or subjected to constraints expressed by a system 
of pseudo -Bo clean inequalities and/or equalities. The optimi- 
zation consists in finding the optimal value of the objective 
function f (x:^ , j and possibly determining one or all 

the optimizing points. 

3.1 Solutions of Linear Pseudo-Boolean Equalities 

Let P-jJ-i + ^ ^2^2 + • • • + ^n^n ^^n~ ^ 

(3.3) 

be the general form of pseudo-Boolean equality where p^, q^ 

(i = 1, n) ands are given constants. can very well 

assume p^ V q^ for all i (otherwise the tem p^^y^^ + q^y^^ will 
be simply a constant and can be brought to right hand side of 
the equation). 


* y^ = (1 - y^)i.e. complement of y^. 
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We will do a little bit of transformation for removing 
complement terms in a non decreasing order. 

Let us set 

■ f if Pi > cii 

X. = < _ (3.4) 

Ui if Pi < li 

then PiWi + = j (Pi“qi) + qi for p^ > q^ 

\ (<li-Pi) + Pi for Pi < cLi 

and bence equation (3.3) becomes 

r^x^ + r^Xg + ... + = t (3.6) 

where r^jr^, ,, ,, r^, t are constants, ri > 0 for each i and 
after reindexing we can suppose that 

2?^ > 2?2 > , . , > > 0 , (3.7) 

Table (3.1 ) discusses 8 mutually exclusive cases concerning 
Eqn. (3.6). 

The solutions of the Eqn, (3.6) which is the canonical 
form of Eqn. (3.3) can be tracked down along the branches of 
tree in I'ig. 3.1 accordingly. Most of the blind alleys 
be avoided by use of Table (3.1). 



TABLE 3,1 





T lie above method will lead to all the solutions of 


the Eqn. (3.6) and from solutions of theEqn. (3.6), the 
solutions of the Eqn, (3.1) can he found. 



PIG. 3.1 
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3.2 Solutions of Linear Pseudo-Boolean Inequality 

' Tlie most general form of a linear pseudo-Boolean 
inequality is either 

p^y^ + q^y -j + P 2^ 2 ^2^2 ^ ^ (5 #8) 

or p-jy^ + q^y-l + ^ ^n^n ^n^n ^ ® (3.9) 

where Pj_> q^, r and s are constants, and we can assume p^ ^ qjl^ 
as usual for all i. If p^, q^ and r are integers, then Eqn, (3.8) 
can he written in the fomn of Eqn. (3.9) hy taking s = r+1. 

In most of the cases p^j qj^ and r are integers, so ,we will 
focus our attention to in Eqn. (3.9). ^he method descrihed 
below for solving Eqn, (3.9) will give solutionsof both 

P^y-^ + q'iyt+. t Pn^n ^n^n ~ ® (3.10) 

p^y^ + q/l + ... + Pn^n + ^n^n > ^ 


(3.11) 



TABLE 3.2 


No. 


Case 


Conclusions 


t < 0 


t > 0 and r-j> • •• 

^ > -t > > .. 


The unic^ue "basic solution is 
= *2 = ... = = ° 
a) J'or every k = 1 , 2 , m: 


n 


= 1 , ■ 

5= .»• 


x, ^ = X, . 
k-1 k+1 


3 t > 0, r.<t (i=1,.,.n) 

n ^ 

and Z r.<t 
i=1 ^ 

4 t>0, r.<t (i=1,...n} 

n^ 

and S r • = t 
i=1 ^ 

5 t^Oy r . ^t (i— i } « • *n) 
n ^ n 
X 3:.>t and Z r <t 

i=1 > 5=2 ^ 


Xj^ = 0 is a "basic solution. 

The other basic solutions (if any) 
are characterized -by the property 

= ••• = 

is a basic solution of 
n 

X r.x > t 
3=m+1 ^ 3 

No solutions 


The uniq.ue basic solution is 


X^ = Zg 


... = x^= 1 


t>0, r . <t (i=1 , 2 , . . .n) 
n ^ 

E r.>t and Z r^>t 
isi1 ^ ■i=^^ 


5=2 


The basic solutions (if ar^r) are 

x,=1, and (x„, ...x ) is a basic 

solution of Z r-X > t - r. 

5=2 3 3 h 

The basic solutions (if any) are 
characterized by the propertys 

either x.=1 and (x„,.,,x ) is a basio 
' n \ ^ n 

solution of .Z 2 ^^x^>t-r^ 
or . ■ / 

x^=0 and (x^, . • • >x„) is a basic 

1 2'n n 

solution of .s r.x. >t, 

3=2 3 3 - V 
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Dgf’lnitioii 3.4 ; 

let S = (y^ , , y^) be a solution of EcLn. (3*9) and 

let I be a set of indices : I 1 , 2, , . , ,n , Let (S,I) be 
set of all vectors (y^, y^) 6 satisfying 

Mr 

7 ^ " ^i i 0 I , 

the other variables y. ( D 0 I) being arbitrary. If all the 
vectors (y^, ,,,, y^) 6 r (S,I) satisfy the inequality (3.9) 
then E (S,l) is said to be a family of solution S of Eq.n, (3*9)» 

If# (I) < n, the family is called nondegenerate* 

We put inequality (3.9) into canonical form as in the 
previous case, 

r^x^ + ^^2 ••• + 2 .'^ (3*12) 

“^n^® (3*13) 

A procedure given below enables to obtain the solutions 
of Eqn. (3.12) grouped into several nondegenerate and pairwise 
disjoint families- of solutions, after this by applying inverse 
transformation we can obtain that of Eqn, (3*9).^ ^^^^^ ^^^^^^^^^^ 

Definition 3.5 : 

A vector satisfying the inequality (3.* 12) is 

called a basic solutibn of (3.12) , if for each index i such 



26 


that Xj_ = 1, the vector (x^ , ... , 0, ..., x^) is 

not a solution of Eqn. (3.12). 

Remark (3.1) ; The solutions of the eqtetLon r^x^ + ...+ = t 

(if any) are also the "basic solutions of Eqn, (3.12), 

The solutions of Eqn. (3.12) are found "by following two steps: 
a) Deteming all the hasic solutions of Eqn, (3.12). 

"b) To each "basic solution associating a certain set of 
indices in such a way that E (S^, 1^^) should "be a 
family of solutions and that the system E . (S|^> I^) 

, should "be ‘complete’ (i.-e., it should include all 
solutions of Eqn. (3.12)X 

Eirst Step : 

a) Determination of the "basic solutions: 

With the help of Ta"ble 3.2, a.11 the "basic solutions ' 
of (3.12) can be found as it was found with the help of Table 3.' 
for Eqn. (3.6). 

Second Step : 

b) Determination of complete system of families of solutions 
of (3 #12): 

To each basic solution S = (x^ > ,. , ,Xj^) we associate 
family of solutions S (S,Jg) as follow s^ Det i^ be the latest 
index for Which x^ = 1, (i.e. x^ = 1 and Xj^ = 0 for all i>i^) 
and let dg be the set of all indices i < iQ. Then S(S, Jg) is 
the set of all vectors (x^, x^) satisfying 
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(xf for i < i 

x' *={ ° ' (3.14) 

i arbitrary for i > in 

V ^ , 

If S^, are all the basic solutions of Eqn. (3.12) 

and if 2. (Sj^j ) (k = 1, ...m) are constructed as above 
k k 

then these families will be ell the solutions of canonical 
inequality (3,12), 

Reverse transformation can be applied to obtain the 
solution of orignal inequality (3.9). 

The solutions of Eqns. (3*15) and (3.16') can be seperatecl 
as 

r^x^ + + r^^x^ > t (3.15) 

r-|X^ + ,,, + r^x^ = t ‘ (5.16) 

For Eqh. (3.16) it is immediate that the set of those solutioi. 
of Eqn. (3.12) which satisfy equality (3.16) will give the 

complete system of solutions of Eqn. (3.16). (See- Remark 3,1 

,/ ' ■ ■ ■ 

Let B. be the set of all basic solutions of Eqn. (3.12)), 
Let be the set of those basic solutions of Eqn. (3,12) 
which are not solutions of Eqn. (3,i6). Let S = (x^, , , , ,x^) 
be an element of B-M’ , Let b be the largest index for which 
= 1, we associate to S , the vectors Rj = • #2-^1 
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for (3 = "b+l, n) defined as follows: 

* 


1 


if i D 


^3i ~ 




1 = if i = 3 


( 3 . 17 ) 


tlie set of all vectors R. for (3 = b+l, ...,n) associated to 

J 

•the different elements of B-M' will Be denoted by M" • let 
M = M‘ U M” . 

Row to find solutions of strict inequality ( 3 . 15 ) follow 
these steps: 

a) Rind M as described above* 

b) Rind the corresponding families of solutions. 

let S(x-] s .. . jX^) denote a pseudo-Boolean equation or inequali 
Definition 3 . 6 : ■ 

The characteristic equation of 2(x^, ,,,,x^) is a 
Boolean equation 

(j) (x^, ...j = 1 (3.18) 

which has the same solutions as i:(x^,,»*x^); the Boolean 
function (j)(x^, ...jX^) will be called characteristic function 
of 2 (x^ » « . * • 



3.5 Determination of Gharacteristic Fimctlon in Various Oase 
a) linear ©.quality.: 


(j) j • * • y ^n^ ~ ^ 


“-t 0^^ 

X ^ X ^ 

1 •••^n 


(3,19) 

2 

wiiere TJ s means that union is extended over all the 

^ -j ? • • *^21 

solutions (a-|,,.,a^) of 2(X|», ,',x^). 


h) linear inequalities: 

*1^ *= 0.J U U Cp 


(3.20) 


where 0^ is the Boolean function oorrespoh- 

ding the to ith family solutions of linear inequality,, 

c) Nonlinear equality or inequality; "■ ■ 

let us consider a nonlinear pseudo -Bo clean 
inequality with unlmowns x.j,.,.Xj^, 

p.|l.^ + , , , + p^I^ > q (3,21) 


where each 1. stands for a certain conjunction (i.e. a 


product of variahles either complement or not): 


• e. 


^k(i) 


‘‘i 


• m ♦ ' 




(3,22) 


*Here “i . „ ' 

Xf - Xj^ If ss 1 
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let us replace the product hy a single bivalent variable 
y^ and solve the resulting linear pseudo-Boolean inag.uality 


Piy-| + 


+ 




> 1 


(3.23) 


where y^, treated as independent variables. 

If ... j^) is the characteristic function of the 

inequality (3.23) obtained as described in (b), then Boolean 
fimction 


e. 


■1 'k(l) 

^ (^'l > ... ^21^ ^ ^ L ’ ... f • • • f > ... 

k(l) “1 



is characteristic function of Eqn. (3.21). Same procedxire can 
be applied for equality, 

3,4 Pseudo-Boolean form of Gharact eristic function 

We have seen that characteristic function of pseudo- 
Boolean system is a Boolean function. However, 1^ 

Basic algorithm for minimizing pseudo-Boolean functions, 
it will be necessary to have a pseudo-Boolean expression of 
chanacteristic function i.e, an expression using only the 
arithmatical operations •*+’» and negation of simple 

variables. These are the following identities for conversion j 
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• • •Ua^ — (a,^+a2+» • •a^) 
n terms 


If (|) t= U O2 UCjjt and if = 0 for all i ^ 3 

tiien 

(j) = + C 2 + ... + • (3.26) 

3.f Minimization of Ifeiconstrained. Nonlinear Pseudo-Boolean 
Itinctions t • 

A vector (x^, £22 is a minimizing. point of 

the pseudo-Boolean function f(X|> if 

f(x^, 1 (3.27) 

H * "if 

for any (x^, ...Xj^) 6B2 , the value f(x.^, ,,,x^) is the 
minimum of f« Brom this definition it follows that 

f(x|, ..., X^) < f (x^5 ...x^, x^^^, ...Xj^ (3.28(1 

for all i =1,2, ...n (3428^ 

Conditions (3.28(1 ))to (3.28(n)) are necessary hut not suffic 


— (a^a2+a-jaj+« »*+a.|a^+a29'^+.*.a23'j^ 

n« terms 
2 

• » • •+a^_^aj^) + ••• • 


n-1 

+ (*"1) ^1^2^3***^n 


n« ■ = It erm 
^n 


(3.25) 
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for (x^ , ,.,,x^) to "be a minimizing point but . they characterize 
local minima of f, 

Basic algorithm ; 

The algorithm is made up of two main stages. In first 
one, the minimum of the given pseudo-Bo clean function 
f(x^, ,.,,x^) is found, while in second all the minimizing 
points are determined. 

Birst stagej let us denote, for the sake of recurrence 

f^(x.j> ***,x^) = f(x^, »..,x^) 
f^ can he expressed as 

f ^ (x-| j • • • = x^g^ (x2> • • • + h-j (x2» ♦ • • 

(3.29(1)) 

Inequality (3.28(1 )) becomes hy using Eqn. (3.29(1)) 

('^■l-5l)gl(2C2» <0 (3.30(1)) 

i !f 

If and Tjr^ are the characteristic functions of the inequality 
< 0 and the equation g^ = 0, x^ can be written as 

■ f ■ V .Jr " / 

x^ = (p^ • * * ’^n"^ “ (^2’***^^n^ ^ P-| f-j (^2’“**^i- 

v(3.3l(l)) 

where p^ is an arbitrary bivalent variable. 
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Tills completes first step. 

At second step 

f2(x2»...,x^) = f^(l}r’(x2s X2,...,Xj^) 

( 3 . 32 ( 2 )) 

Since f^ is a pseudo-Bool ean function, TjT^ should be in pseudo- 
Boolean form,' Hence 

f2(x2, .,.,Xj^) = (X2J . g^(x2» + li^(x2»...,y 

Continuing this way, we obtain at nth step a function f^(xjj) 
which is written in the form 

+ ‘"n (5.29(n)) 

where g^^ and h^ are constants. 

We have to solve inequality 

- ^n) Sn i ° 

The solution is 

X ss ( p ) 

r 1 if < 0 

Where = J O' if > 0 

The ■ miniittum of orignal function f is 
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is given by 

if gn < 0 
if Sn > 0 
to an end. 

Second stage: Tbe various minimizing points can be 
determined by recursion as follows: 

^n ~ ^^n^ 

=^n-1 = fn-l 
♦ 

(Pf? • • • > ’ ^n^ 

Recently many new methods for minimising the pseudo 
Boolean functions have been suggest 


fv.- = f ^ 
min -^m-l 


where the -constant t 


n 


^n 




0 


!fhe first stage comes 



CHAPTER IV 


PSEUIX)-B00LEAH EORltHLATIOH 

This chapter is devoted to the development of the 
methods for getting the optimal results for all the three 
problems. This has been done by reducing the task to that 
of minimization of unconstrained nonlinear pseudo -Boolean 
functions. Proofs have also been given to ensure that mini- 
mization of these unconstrained nonlinear pseudo-Boolean 
functions will give the optimal results. 


4.1 The m-Vertex Median of a G-raph; 

let us recall of our notations first: 


G : Given grapth with a veirticesj 

J ftLstance natri^c ’(u x a), represents the minimum 

distance between vertex v. and v. . 

X J 

X : nth order diagonal matrix with vertex wei^ts on 

the diagonal. 

H : Wei^ted distance matrix, h^^^ representing 

minimum weighted distance of vertex v_. from v^. 


H = B X = 1 I = I Xj I 

We have to find a source subset ( ^ ~ m) such that 
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1^* < lip 


for p6 (jjj) 


where n 

h = 2 

p 3=1 v^evP 


Min (h. ^ ) 


(4.1) 


Define the following variables; 


74 = 


1, if a source is located at vertex v^ 


= 0, otherwise 

t .. = 1, if there is a source at v. and v . served by v. < 

-L. J J. 

= 0, otherwise. 


This problem (say P) can be written as follows; 


Problem 

P 


n n 

Min 2 2 t-.^ ii.V 

i=i 

n 

2 y . = m 
i=1 


n 

2 t. 


i=1 


ID 


(for all j=1,,,*,n) 


and y. = O mr^k-t, . - 0 for all j = 1,2, n 

1 ' ID 


(4.2) 

(4.5) 

(4.4) 

(4.5) 


One can easily makre out that t . . are dependent variables 
as for a given Y vector, of yj_*s, the corresponding T 
matrices of t,-4; ’s can be easily determined. 
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Let us define 




0 

0 

1 

1 


if h. . < li, . 

13 . ko 

if h. . = In . and i < k 

13 . iCD “ 

if h. . = K. and i > k 

13 -*^3 

^id > ^kd . 


(4*6) 


Express t . . as follows: 
13 


n 


^id 


= ^1 ■ Vi)c3> 


(4*7) 


Let us consider tlie following problenj 


n n 


n 


Problem 

p f , 


(4.8) 

(4.9) 


n 

E y . = n 
i=1 ^ 


Lbeorem 4.1 1 fbe minimum of problem P is equal' to tbe 

')f ' ' ' 

minimum of Problem P*. If (y-| , , y^) is an optimal solution 

* • 

to problem P* and if tj_^ the corresponding ■values of t^^. ^ 

(by expression 4.7 ) » (y-j > . • ^ "^11 * • • • >^nn^ optimal 

solution to P. 
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Ihe proof of tlie theorem will "be followed hy the 
proofs of the following lemmas s 

Lemma I t If S 5 = is a feasible solution of 

the Problem P* then the corresponding solution S* = (y^»***>yj^j 
t^^j t^) will be a feasible solution to Problem P 

(tji^^ *s are determined according to expression(4*7))* ’ 

Proofs We will prove that S* satisfies all constraints of 
the Problem P thus making it feasible, 

S’ will satisfy constraint (4*3) as S satisfies 
constraint (4*9)* 


t^^ in S* are expressed as follows s 


= yi - yk "'ikP 


In the above expression if y^ = 0 , then t^^. *= 0 
for 3 = 1,,,,.,n. Therefore, solution S* satisfies the cohs- 
traint. set ,(4*5) too'*! 

Now %e have to show that S’ satisfies the constraint 
set (4*4)' too. 


Let us define Ij^ *= | ^ 1 ^ 


Prom above expression for t^^ and definition of 

according to ( 4 * 6 ) we will prove below that for eada 3 , there 

exists one and only one s^ such that ^ j *= 1 and f or that 
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» "tile following are tmie: 


yg. = 1 


and 

eitlier 




for (4.10) 


or 




We have to, first show that if systeir. (4.10) ie satisfied 


then t _ . =1 


n 

■to 4 =70 ^ Cl - 7]. V Jj. 4 ) 



V 


as either 
or 


Hence t 



2S 


H 





s.lcj “ ° 

.:s 1 



For any 4, there always exists one s^ satisf7"ing syst^ 
(4.10). 



4-0 


We will prove that for every 


= 0 -for p ^ s^. 


If P0I, then y = 0 and subsequently t . = 0 . So we have 


to prove 


t j = 0 for and p V Sj 


jFrom system (4*10) 


either hj,j > h^.j 




(4.11) 


System (4.11) says that v^ = 1 


t . = % (1 - v^ ^.) for all 

1< , 

peij^ and p je s^ • 

== 0 

(as for k = s . , ’^pe^ = ^ 


Therefore S t_j = S t ^ + S t j + t- ^ 

p^t P0I^ peij, ^3 s^3 

I^s. 


(4.12) 


’s4#3 


1 
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Now Eqn. (4f12) .is true for every 3 , therefore it shows that 
S* satisfies the constraint set C4»4) too* 

Hence S* satisfies all the constraints of Problem P, 
thus proving lemma . 4 ,' 1 * 

let (t^^, y^, y^) be any feasible 

solution to Problem P and let f* denote the value of the 
objective function of P at this point. 


let us put 



* ^ 

= y^ Ci 

^ 1&=1 

* 

ij , 

" ^ija 


^ * # 

• • • > 

11’ 

• ♦ ♦ f ^ *1 ^ 





function of P at this point. 


Iiemma 4,2 ; f C f 

* n n * 

Proof: f = 2 2 h.. .• t . . 

3=1 i=1 


n n n 

2 h j , + 2 £ (h» j 

3=1 3=1 i=i ^3 




n 

3*1 3 


Hence leEma 4*2 is true. 
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It follows directly from lemma 4.1 and lemma 4.2 
that the minimmn of problem P* is also the minimum of 
problem P. 

Thus we have reduced P to P'. 
n 

let Y = E Max (M .) 

3=1 i=1,2, 

• • «i « 5 ^ 

lemma 4.3 : Any solution to the problem P* or P will never 
give the value of the objective function more than Y, 

Proof:- Obivious, 


Consider problem P ” 

n n n ■ n 

Problem Min E E 7 ^ % (1-y-f^ v.^ . ) h. . + Cm+1 ) ( 1 

D=1 1=1 1&=1 ^ ^ i=1 ^ 

. p tf 

(4.13) 

Theorem 4.2 : The minimum of the problem P’* will also be 
the minimum of problem P’, 

Proof of this theorem will follow the proof of 
following lemma; 

lemma 4.4: Any optimal solution S which gives minimum value 

• " n 

of (4.13) (problem P“) satisfies S 7 ^ - m, 

i=1 

n ■ 

Proof: If 2 y. T^ m, then value of (4.13) will be greater 

than Y contradicting the fact that S gives minimum value 

n V • 

of (4. 13 ) I because any solution satisfsring . E y. = m in (4.13) 

i=1 
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will always give the value of (4.13) not exceeding Y, 

Hence lemma is true, 

The proof of the theorem directly follows from 

n 

lemma 4.4j "because if S Y. = m then (4.13) and (4.8) 

i=1 ^ 

becomes same; satisfying constraint (4.9) too, Therefore 
the minimum of will be the minimum of P'-too. 

It follows directly from Theorems (4.1) and (4.2) 
that the minimum of P‘*will be the minimum of P too. Thus 
we have reduced the problem to minimization of an unconstrain 
nonllnaor pseudo-Boole an function. 

Hammer has proposed Bessie algorithm to minimize 
these unconstrained, nonlinear pseudo-Boolean functions. 

In Chapter Y, some important points are mentioned to 
code this algorithm efficiently on computer. 

Some other ideas for minimizing these types of func- 
tions have also been proposed recently^^^^\ 

formulation and proofs for the third problem will 
help in better understanding of formulation and proofs for 
the second problem, therefore the third problem follows first 

4*2 finding Minimum Humber of Bn erg ency Service 
let us look again at the problem, we have 

s. 1 the maximum weighted distance alio?; ed that seperates 
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Vertex from its nearest source vertex. 

We have to find the minimum number of sources 
required to meet these distance standards with other things 
as usual. 

If we define p • • as 

J 

= 1 If h. ^ < 

-LJ 

then this problem can be expressed mathematically as follows 



: n . 

Problem Min E y. 

Q , 

n 


( 4 . 14 ) 


Wi > 1 

for (3 =1,2, n) (4.15) 

Consider the following problem . 

n n n ■ 

Problem Min 2 y. + (n+l) E 2 (p<h T.) (4.16) 

qV i=1 ^ D=1 i=1 

theorem 4.3 1 Ihe minimum for the problem is also the 
minimimi for problem Q, 


2 P 
i=1 


Id 
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T3ae proof of tiie above theorem will follow from the 
following lemmas: 

lemma 4 > 5 : The minima for both Q and will never be 
greater than n. 

Proof: T airing solution S with all y^s equal to one, we 
observe the following: 


(i) Solution S is feasible to problem Q, as the 

constraint set (4»15) is satisfied, as — t 
for (o’ = 1 •>»■)• 

(ii) The value of the objective function of Q 
corresponding to solution S is n, 

(iii) Solution S is feasible to QS as Q‘ does not 
have any constraint, 

(iv) The value of objective function of Q* correspondir 
to solution S is n. 


The proof of the lemma follows directly from (i)» 
(ii), (iii), (iv). 

Lemma +.6: -toy optimal solution S of problem Q* satisfies 



(4.17) 


Proof! If system (4.17) is uot satisfied for any of j then 
the objective function (4.16) will be greater than n and 
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hence contradicting the fact that S gives mininanm of Q*» 

Lemma 4 » 7 t Any optimal solution S of the problem Q‘ will 
satisfy the constraint set ( 4 «’ 15 ) i.e. 


n 


2 P.. y. > 1 
L=1 -3 -L 


for (j = 1,2, .,n) 


Proof; By lemma 4.6, S will satisfy 


n 


S (p„ y^) = 0 


i=1 


for (o = 1 , . . . >»■) 


This implies for every 3, there should exist atleast 
one i , such that 


n 




K Hi ^ ^ 


1=1 


n 


Therefore satisfies ^ 2 ^ Pij 

for (3 = 1,2, ...jn) 


Lemma 4 . 8 ; Any solution, satisfying the constraint set 

(4.45 ) i.e* 


z 5i3 n i 
1=1 


f or ( 3 = 1 , 2 , » • • , h ) , 
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n 

satisfies % (p. . yT) = o 

foi'(3~”l>2j*»*jP3 

Proof: Obivious by observation. 

It follows from the above lemma that for all feasible 
solutions to problem Q, the problem is 

n 

Min £ y. 
i=1 ^ . 

and lemma 4.6 says minimum of Q' will always satisfy 
n n 

2 ^ CpJTl^) = 0 , 

3=1 i=i J 

Prom the above two statements and lemma 4.7 V i't follows , 
directly that the minimum of Q* is also the minimum of Q, 

Thus once again we have reduced this problem to the 
minimization of an unconstrained, nonlinear pseudo-Bool ean 
function. 

The second problem follows: 

4.3 The Modified m~?ert ex Median of a Graph ^T^^ 

As already mentioned in 2,2, in this modified m-yert ex 
median problem we have to satisfy the distance constraints 
of the third problem too, when solving the first problem, 
m-vertex median. 
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The problem can be expressed mathematically as 
follows : 



n 21 

n 


Problem 

R 

Min 2 ^7. 

0=1 i=1 ^ 


h^^ + (Y-fl) 



n • 2 




( ^ y. - m ) 
i=1 ^ 

(4.18) 

such that 

n 




I y. > 1 

i=1 ^ - 




for 

(o — 1»2,,», ,n) 

(4.19) 


It may happen that there exists no feasible solution 
to modified m-vert ex median of the given graph &, This 
depends on s.’s and m for a given H for the graph G-. Kow 
even if there does not exist a feasible solution to mr-vertex 
median of &, there will exist a feasible solution to R, but 
minimum of R will be larger than Y. We already know minimum 
of R Y^ill never be larger than Y if there exists a feasible 
solution to modified m-vert ex median. 

Hence, this information should be kept in mind? if 
minimum of R happens to. be more than Y, then there does not 
exist a feasible solution to the modified m-vert ex median. 

Consider problem R’ as follows 
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Problem 

Rf 


n 

Min s 
D=1 


n 


Z 7. 
i=1 ^ 




+ (V+1) 


n 2 ^ ^ 

( E y. - m) + (v+1) £ 2 

i=1 ^ 3=1 i=1 


^yi I’i j '> 


Theorem 4«3 : If there exists a feasible solution to the 
modified m-vertex median of the given graph G-, then the 
minimum of R’ Will be the minimum of R* 

The proof will follow the proof of the lemma given 


below , 


Lemma 4*9 ? If there exists a feasible solution to the modi- 
fied m-vertex median of the given graph G, then a solution 
S giving minimum of R’ will satisfy the constraint set 
(4.19) of R, 


Proofs Row if there exists a feasible solution to modified 
m-vertex median of G, then the minimum of R’ will not be 
greater than V, hence S will satisfy 


n 

71 


i=1 


(Jl Pij) = 0 


for (3 


1 , ...,n) 


system (4*20) shows 
^ ^i^i3 

^ for (3' =1, n) 


(4.20) 
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Hence S satisfies constraint set (4»19) of R i.e, 
it is a feasible solution of R* 

This proves lemma , 4 , 9 . 

Lemma 4«10: Any feasible solution S of R will satisfy 


Since 


Hence 

Ibis proves lemma 4.10. 

The proof of theorem 4.5 follows directly from proof 
of lemmas 4.9 and 4 • 10. 

Here also we have reduced the problem to minimization 
of unconstrained, nonlinear pseudo-Boolean function. 

While we minimiz;e these unconstrained, nonlinear 
pseudo-Boolean functions we-get not only single optimal 
solution for the problem, but also all the optimal solutions 
existing for that problem. 


n n 


.2 n: Cy, P..) = 

1=1 i=-l 1 


0 


n 


2 y. Pii > 1 

i=-| ^ 


for (3 — 1 , . . . ,n) , 


n 

% 


1=1 (n = 0 


for (0 = 1,..,,n) 


n n 


2 % (y. P. .) = 0 

j=1 i=i ^ ^3 



CHAPTER V 

COBfIPUTATIOHAL ATTRIBUTES 

In this chapter principles of computer code for 
Basic algorithm to f ind aH m~vertex medians of a graph are 
described. Also the computational attributes of various 
heuristic algorithms for the above problem are compared. 
Summary of the results obtained on IBM/7044 for problems 
of various sizes and the conclusions made for them are 
also given. 

5.1 Coding of Basic Algorithm 

In Chapter III? Basic algorithm to minimize uncons- 
trained peeudo-Boolean functions was described. In. brief, 
it contained essentially the following: 

let f (y.|> ... ?yj^) be the function to be minimized. 
Consider the following system 

Cs.id),) 

for (i 1,2,...,n) • . , vx 

(5.l(n» 

If we successively solve the inequalities ( 5 . 1 (l)) to (5,1 (n) 
introducing in each inequality the parametric solution of 
the preceding one, we will obtain all the minimizing points. 
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Jliere are . essentially two stages. In first stage we find, 
ninimuin of f and in the second stage all the minimizing 
points of f , 

First stages . For the sake of recurrence we write 

f -i (y-j j y^^) = f (y..^ 5 • » * jy^^) (5*2) 

and express f^(y.,,...,y^) = y.jg^ (y 2 > • ♦ • jy^^) + h^ (y 2 » •• ‘^y^) 

(5.5) 

. f ' ft 

We find and fi the characteristic functions of g.^ < 0 
and g^ == 0 and then we get f 2 as follows j 

f2(y2’ “ ^1- ^^1 ^^2’ •• ^2^****^ri^ 

I 

By converting in pseudo-Boolean form 

^2^^2’ ‘ **’^n^ “ *’^n^ g-j (y2» * * * + h^ (y2? * • • 

We store the value of y.^ in terms of (y 2 f...,yj^) given hy 
the following expression (5.5)> "'''^kich will he used in the 
second stage to determine all the minimizing points of f. 

y^ =!= 'f.jCpif y2> • • • fy^i) (y2» • * U p.^ (y2» • >ryin 

Where is an arbitrary Mvalent variable. Ihis completes^ 
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function fg in. the same way 
after nth. step we have 

C5.6) 

C5.7) 

where is an arbitrary bivalent variable, 

The first stage comes to an end. 

Second stage: Introducing the values of y^ (obtained above) 
in the expression for we obtain the corresponding values 

of and then introducing values of y^^ and in the 

expression for yj^_£ we get the values of yj^_2* Continuing 
in this way we obtain all the minimizing points of f. 

Considering the problem of m-vertex median, we have 
to minimize the following unconstrained, nonlinear pseudo- 
Boolean function, 

, . . . ,y^) = l ^.2 Yi J Cl-Ji, ) Hij + (T-H ) ' 

(■5.8> 



first step. 

We proceed now with the 
as with f Continuing this way 

^nin = f = f ClT^) 
n+1 n '•^n^ 

and 

’'n = (Pn> 
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ilie nimiber of toms involved and number of elements 
in tlie terms for expression (5.8) are highly undeterministic, 
depending upon the problem. Using 'SLIP»* routines for 
coding Basic algorithm will be an -elegant way to minimize 
above function. 

One way to solve this problem is to generate the 
whole expression (5«8) from given H matrix, store it in the 
core memory and then apply Basic algorithm as stated above. 

In this way the memory requirement will be too large limiting 
the size of problem which can be solved on a digital computer. 

We observed in the description of Basic algorithm that 

at any step i in the first stage we require only g^ for 

f n , 

getting and This observation suggests an elegant 

way of utilization of memory by not putting the whole expr- 
ession (5.8) in the core to start with, but generating it 

in parts as and when the need arises. this procediire 

■¥: 

we will be minimizing the following function f : 

f*= f - (v+1) m^ (5*4) 

2 ' 

Since (V+1) m is a constant termj the minimizing points 

■ ' # . 

of f will be same as that of f and 


* f . 2 

*= + (7+1) m 


( 5 . 9 ) 


routines for list processing are available on IBM/7041 
at I.I.T, Kanpur, India. 



55 


At first step we will genera.te all those tenas whioh. 
have y-| as an element, let us denote if as f ^ . f^ is 
given by the following expression' 

t n n p ja 

J0¥* ' . 


n n 


3^1 J 2 '*-3 "33 Z 2 ^^^7^ ■ ^5.10) 


We will make f^ = f^ and find g^ by Eqn. (5*3). 

, ,t . Tt 

Subsequently we find and ijT from this g^ as usual, 

Wg find £ 2 ^ ^^2^ follows* 


ir f 

^2 ^y2>***»yn^ ~ » (3^2’ * • • '^n^ • • • *^n^ 


(5.11) 


Now we generate f^ by, following expression * 


yj = ^2,1 3'23 ’"2k3^ 


- ^2 J3 hj h "123 " ( 3 -T„ 


+ (V+1) 72 


[0- 


2m) + 2 E y 
s=s3 ® 


(5.12) 


* The expression (5.11) can be obtained by putting r's= 1 
in expression (5*13). 

** The expression (5.12) can be obtained by putting r =2 

in expression (5.13). 
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let us construct the function f 2 as follows: 

t IT 

We proceed now v/ith f 2 in the sane way as with 
f^. At step r we will find 


f IT 

fr = + fr 


where 




y 2 h . 7t (1-yj^ 
^ ; 5;=1 ^3 




n n n 

- y 2 2 h.. y. 

,1=1 i_r+1 ^ k^r+i 


+ (y+1 ) y:p I (l-2ni) + 22 y 1 
n_ s=r+1 J 

( 5 . 15 ) 


Continuing in this way after nth step we have 


^min “ ^n+1 ^n+1 % V 


( 5 . 14 ) 


yn ^ % 


( 5 . 15 ) 
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v^here is an arbitrary bivalent variable. 


^min ” ^min ^ 


After this, all the minimizing points of f will be deter- 


mined aa usual. 


These observa-tions direct' us to develop an algorithm 


as follows: 


Step 1: Generate f^ by e3q)ression (5,10) and make 

. t . ti 

Step 2: Find g^ and from that corresponding > f-j » f| 

« . 

and f 2 . Save 

Step 3i t*--2 

Step 4; Generate f^ by expression (5* 15). 

Step 5: f^<— f^ + f^ 

Step ^5: Find corresponding g^ from fj.* ^nd from g^’ the 
corresponding j ^r+i • rr* 

St ep 7: If r < n then r«— r+1 and go to Step 4. If not 
go to Step 8. 

* tf 

Step 8: Vl 

step 9: , 
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Step 10j Dot ermine all the minimizing points from, 

1^1’ f2^****^n described in Stage 2 of the 
Basic algorithm. 

At Step 9 > "We get the minimon value of objective 
function i.e, and at Step 10, we get all the minimizing 

point of f, 

A similar approach can be adopted for the second 
problem i,e* modified m-vertex median and the third problem 
i.e. finding minimum number of emergency services. 

A code for the algorithm described above was pre- 
pared using *SIiIP^ routines available on Oonputer IEM/7044* 
This code has mainly three routines whose functions are 
given below s 

Ci) Routine OPTAL ?- This routine creates list structures 
IPRO, ISTO, lYAR and ISVAR. The lists ISTO AND IPRO are 
used to store f^ in parts at any step r. list IPRO contains 
expression for g^ (one sublist of IPRO representing one 
term of g^) before it calls the Routine BALGO, G?he Routine 
BAIGO and subseq.uently Routine AIiSOIN do necessary things 

t ' ' ' ** 

to store t and in lists ISVAR and IV AR respectively, 
and also update IPRO, Now we increase index r to rH-1 and : 

repeat the whole procedure. When value of r becomes n this 
Routine OP TAD finds f^^^ , i.e, the minimum value of function 
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f , and all values of variable y^. All the minimizing points 
of f are now found by processing lists IS7AR and IVAR* 

Bp ut ine B AIiGO ; It processes list IPRO to get all the 
solutions of ^ < 0 and ^ = G, and stores them in list 
ISOM, Out of all sublists of ISOM, those sublists which 
have a mark correspond to the solutions of g^ =0, 

After this it calls Routine AIiSOM, Finally it 
updates IPRO and returns back, 

r 1! 

(iii) Routine ALSOM ; Ihis routine finds and from 

list ISOM and stores them in lists ISVAR and ITAR respective .y* 

. t 

is converted into its pseudo-Boolean form and put in 
list HEMP, Now it multiplies lists IPRO and HEMP and 
puts product in list IMHMP, After this it returns back. 

The flow-charts for these routines are given on 


the following 3 pages,. 



Routine OPlEili 
f " ST ART > 


CALL LIST 
CALL LIST 
CAIL LIST 
CALL LIST 

r = 1 , f 


[IPRO) 

ISTO) 

IVAR) 

,ISVAR) 


Scan IPRO, keep only tliose sub- 
lists in IPRO whicb. Rave y^ as 

an elenent, RemoTe other sub- 
lists in IPRO and store then in 
ISTO 








Eoutine BilGO 


STikET 


Find all solutions of fron IBRO 

and store then as the sublists of 
I SO IN, Mark those sublists of ISOIN, 
which correspond to solution of g„ = 0, 


STiiRT SCiiMING BIST ISOIN 


Advance by one sublist 


Tes ■ 


Header of list ISOIN \ ^ Ho 
encountered? / 











Routine ALSOIN 
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5*2 Comparison of Gommtatlonal Attributes 

Both Maranz ana* s algorithm and One optimal 
algorithm do not ensure optimal results. It is quite logical 
to think* what will happen if we apply these algorithms 
alternatively as follows: 

Let he any- solution chosen at random. With 
as an initial solution let and 0^ he the final solutions 
hy Maranzana's algorithm and (One optimal algorithm respective!'" 
Algorithm ; 

Step 1; 

Step 2> Apply Mar anzana*s algorithm taking A as initial 
solution to get B as final solution. 

Step 5s If B is same as A, go to Step 6, If not, A **"~B 
follow Step 4* 

Step 4; Apply One optimal .algorithm taking A as an initial 
solution to get B as final solution. 

Step 5; If A is same as , B, go to Step 6, If not, A-^-B, 
go to Step 2, \ 

S' ■ ' ' 

Step 6: Stop, B will he the final solution hy this 

procedure* 

Among the several conputer runs for comparing 
computational attributes, in most of the cases it was found 
that B and 0^ were same. 
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In the above algorithm if we interchange the 
position of Steps 2 and 4» then in all computer rims B was 
found to be same as G^, 

Table 5*1 summarizes the results of various computer 
runs for comparing the computational attributes of Maranzana^s 
algorithm and One optimal algorithm, Bor each entry in 
columns 3 and 4 mean times ore taken over ten different random 
initial solutions. 

If for each entry 

M : Mean value of the ob;} active fimction by Maranzana*s 

algorithm taken over ten random initial solutions 

R-j j * **, 

0 j Mean value of objective function by One optimal 

algorithm taken over ten random initial solutions 

R-j > • * * 

Maximum value of objective function by Maranzana'*s 
algorithm taken over ten random initial solutions 

R-j 5 ♦ • • i R j Q • 

0 : Maximum value of objective function by 'One optimal 

X ' \ ^ 

algorithm taken over 10 tandom initial solutions 

R-j » « • • jRjO* 




65 


Then in colurm 5 

Mean error % = x 100. 

In colinnn 6 

Max error 7* 

Oonclusions : 

(1) The final solution, by Maranzana’s algorithm was 
found to be highly sensitive to the initial solution, 
while by <kie optimal algorithm in most of the cases 
the final solution was found to be constant* 

(2) Pot a given initial solution, the final solution 
obtained by One optimal algorithm was never found to 

be inferior to the final solution obtained by 

Maranzana^s algorithm^ 

(3) Any One optimal solution taken as initial solution 
for Maranzana's algarithm was found to give the same 

final solution. 

(4) Since mean error (varying from 67. to 687*) and 
maximum error (varying from 227 * to 206*/,) were varying 
between very wide ra^es it is very rislsy to solve 

the problem by Maranzaha*s algorithm, althou^ time 
taken by it may be lesser than that by One optimal 
algorithm. 


M - 0. 


X 


X 100 


X 



No, of verti- 
oes 
n 


No, of 
medians 
n 


Mean comp— Mean compu— 

utation ti- fation time 

me for Mara** for One— opt, 
nzana's algo— algorithm, 
rithm. Time Time unit 
miit 1/60 sec. 1/60 .sec. 


Mean 

error 

V. 


Max, 

error 


10 

10 


2 

4 


5 

10 


8 


37 


10 

5 

13 

15 

3 

13 

15 

6 

24 

20 

2 

13 

20 

4 

27 

20 

6 

40 

20 

8 

52 

20 

10 

61 

25 

2 

22 

25 

5 

60 

25 

10 

107 

30 

3 

44 

30 

6 

108 

30 

9 

172 

30 

12 

186 

30 

15 

240- 


15 

49 

112 

24 

84 

194 

38 

26 

105 

109 

41 

88 

• 44 

10 

24 

166 

12 

22 

287 

12 

29 

562 

26 

74 

546 

38 

57 

69 

9 

• 26 

387 

13 

49 

1024 

44 

104 

292 

12 

63 

1015 

33 

92 

1785 

49 

75 

2538 

68 

206 

2560 

67 

143 




COMPUTATION TIME ( UNIT - Vgo SEC- ) 



NO. OF VERTICES = 10 
X represent time taken 
maranzana's method 

O REPRESENT TIME TAKEN 
ONE OPTIMAL METHOD 



BY 

BY 


COMPUTATION TIME ( UNIT-Ugg SEC ) 

bi? ^ -C^ cn 
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